Density fluctuations of a hard-core Bose gas in a one-dimensional lattice 
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For a hard-core Bose gas on a one-dimensional lattice we find characteristic oscillations in the 
density-density correlation function. Their wavelength diverges as the system undergoes a continu- 
ous transition from an incommensurate to a Mott insulating phase. The associated static structure 
factor vanishes as the Mott insulating phase is approached. The qualitative picture is unchanged 
when a weak confinig potential is applied to the system. 
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Recent experiments on one-dimensional (ID) Bose sys- 
tems have opened an exciting new field of physics. 
This will provide us with an extended and deeper under- 
standing of the special properties of strongly interacting 
particles at low dimensions. From the theoretical point 
of view ID systems are easier to treat in comparison with 
two- and three-dimensional systems but also prevent us 
from using conventional mean-field methods. 

A well-known fact of one-dimensional many-body 
physics is that a continuous hard-core Bose gas is equiv- 
alent to a free Fermi gas 0, [f| . The role of the hard-core 
interaction of bosons is played by the statistical prop- 
erties of the fermions, related to Pauli's principle. This 
idea can be conveyed to ID lattice systems. 

In the following we will adopt an approach to the statis- 
tics of directed polymers [6| in two dimensions. The 
analogy with this classical statistical problem is based on 
the observation that the world lines of a grand-canonical 
ensemble of bosons are equivalent to directed polymers, 
random walks or fluctuating flux lines 0, El • This can be 
formally expressed by the fact that the partition func- 
tion Z of the grand-canonical ensemble of these systems 
is identical. 

For directed polymers in two dimensions it was shown 
that Z can be expressed as a determinant @. Thus the 
partition function of hard-core bosons in d = 1 reads 



Z = dct R 



(1) 



where R is diagonal with respect to the Matsubara fre- 
quency u> 



f-l + e i^f} ai 



e™T)a 2 . (2) 



T is the shift operator along the ID lattice (Tf(r) 



1)), and the <jj are the Pauli matrices 
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This model describes the tunneling of bosons with rate 
J > between nearest neighbors, expressed by the shift 



operator T and its inverse f . J is measured in units 
of the lattice potential. The 2x2 structure arises from 
the fact that particle exchange between neighboring sites 
through simultaneous tunneling is prohibited. This re- 
duces the translational symmetry to sublattices with ev- 
ery second site, where r is the position on a sublattice. 
C > is the fugacity that controls the density of bosons 
in the system. The fugacity is not directly accessible in 
the experiment but only indirectly through the density 
J). Therefore, physical quantities should be mea- 
sured as functions of n and J. An additional potential, 
superimposed on the optical lattice, is described by a 
space-dependent fugacity ( r . 

Physical quantities can be derived from the matrix R. 
For instance, if we are interested in properties at zero 
temperature we can use the integral with respect to the 
Matsubara frequency 



2^ 



f-ZT> 

Gf f j*f — / )t".t' 

to evaluate the local density of bosons as 

7l r = 1 + £ r ^G r ^ r 



(3) 



(4) 



It should be noticed that G T y is not the Green's func- 
tion for the propagation of an individual boson in the 
system but a correlation function between two bosons. 
Therefore, it is not possible to evaluate the momentum 
distribution of bosons 0, S IE S E3, 113, E3, E3, E3 from 
this expression. On the other hand, the evaluation of 
the density and the density-density correlation function 
becomes a simple task with the matrix G r y [|J . 

Another interesting quantity is the correlation function 
of the density fluctuations 



C r y = -((n r - (n r )) (n r > - (n r /))) 

(CrCr' ) G r ,r' Gr' ,r 



(5) 



from which the static structure factor as a function of the 
momentum k can be obtained by Fourier transformation: 



—ikr 



(6) 
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FIG. 1: The zero-temperature phase diagram of the model 
shows three phases: a Mott insulator (MI), an incommensu- 
rate (ICP) and an empty phase. For J > Jo = 1/2 there is 
no MI. 



N is the number of sites of the ID lattice. 

The local density n r and the correlation of the density 
fluctuations can be directly measured in an experiment 
[l5| . This motivates the following study of these quanti- 
ties for a translational-invariant system as well as for a 
system with a weak parabolic potential. We will compare 
the results in the incommensurate regime near the transi- 
tion to the Mott insulator and discuss their characteristic 
properties. 

(i) Translational-invariant case. For constant fugacity 
Q r = C the system has translational symmetry on the 
sublattices, and G r y can be calculated analytically Q. 

Figure ^ shows the zero-temperature phase diagram 
of the model. The particle density is constant in space. 
Three phases can be identified: an empty phase with n = 
for C < 1/(1 + 2J), a Mott insulator (MI) with n = 1 
for £ > 1/(1 — 2 J) and J < 1/2, and an incommensurate 
phase (ICP). For J > J = 1/2 the system exhibits no 
MI phase. The density in the ICP can be calculated from 
Eq. (@J and gives 



n = 1 



1 

27 



k =f (k* - tt) 



(7) 



spectively, and k, k* are given by 

(i + r 1 ) 2 



k 
k* 



arccos 1 



2CC- 1 + J 2 ) 
(i-C 1 ) 2 



2J 2 



- 1 



(8) 
(9) 



The transition from the intermediate to the Mott insu- 
lating phase at the critical fugacity £ c = 1/(1 — 2 J) is 
continuous. 

To investigate the behavior of the system near the Mott 
transition we calculate the correlations of density fluctu- 
ations asymptotically from Eq. (J5J as 



sin (k*r) 



(10) 



After a Fourier transformation we obtain the static struc- 
ture factor as 



S k 



2^ ( k * h k ) 

2^2 (k* + \k) 




< k < 2k* 
> -k > -2k* 
\k\ > 2k* 



(11) 



where =p correspond to the cases £ > 1 and £ < 1, re- 



These quantities are shown in Fig. |2 a,b for two values 
of the tunneling rate J in the ICP. They vanish as the 
MI phase is reached due to the fact that the MI exhibits 
no density fluctuations. The correlation function of the 
density fluctuations shows significant oscillations in the 
ICP. Their wavelength A = 2ir/k* determines the charac- 
teristic length scale for density fluctuations and diverges 
as the Mott transition is approached. S k is nonzero in 
the interval [—2k*, 2k*} and falls off linearly to both sides 
from k = 0. The slope of the static structure factor does 
not depend on the tunneling rate. 

The dependence of k* on the density for varying tun- 
neling rates is depicted in Fig. At low densities k* 
increases with increasing n until it reaches its maximal 
value of 7r at a certain density, where it shows a cusp for 
J < oo. For J < Jo the system can undergo a Mott tran- 
sition, where k* vanishes at n = 1. Otherwise it stays 
nonzero. 

Note that these features of hard-core bosons in ID dif- 
fer fundamentally from those of an ideal Bose gas in ID. 
There Co,r decays exponentially and its oscillations show 
no characteristic behavior. 

(ii) Parabolic background potential. A parabolic po- 
tential can be expressed as a spatially varying fugacity 
CT 1 = + fir 2 , where ft determines the strength of 
the potential. In this case G r y cannot be evaluated sim- 
ply by a Fourier transformation, since the translational 
invariance is broken. We have calculated the local par- 
ticle density the correlations of the density fluctuations 
and the static structure factor by inverting the matrix R 
numerically on a lattice with N = 500 sites. 
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FIG. 2: Correlation function of density fluctuations Co.r and static structure factor Sk for different tunneling rates J. First 
row: constant background (£ _1 = 0.3, f2 = 0). Second row: parabolic background (C" 1 = 0.3, £7 = 3- 10 ). Tunneling rates: 
Ji = 0.3506 (solid lines) and J2 = 0.3504 (dashed lines) 




FIG. 3: Characteristic wave vector k* as a function of the 
density in the translational-invariant case. Curves are plotted 
for different values of the tunneling rate: J = 0.2 < Jo (solid), 
J = 0.5 = Jo (dashed), J = 0.8 > J (dotted). 



FIG. 4: Local particle density for parabolic background po- 
tential (C _1 = 0.3, fi = 3 • 10" 5 ). Development of a Mott 
plateau in the center of the trap (r = 0) as the tunneling rate 
J is decreased below Jp w 0.35. 



The local particle density n r is shown in Fig. 0] for 
different values of the tunneling rate J. The density is 
symmetric around the minimum of the parabolic poten- 
tial at r — with a maximum at the center (lfjl |. It is 
suppressed as the potential becomes larger with increas- 
ing distance from the center of the trap. As the tunneling 
rate is decreased the distribution of the particles along 
the lattice becomes narrower and the density is shifted 



upwards. When J reaches some value Jp, we observe 
a region with local particle density n r — 1 developing 
symmetrically around r = 0. 

To investigate the development of this plateau we have 
evaluated the correlations of the density fluctuations 
Co, r together with the associated static structure fac- 
tor. These quantities are depicted in Fig. [2]c,d for two 



4 



values of J > Jp. The correlation function of the den- 
sity fluctuations exhibits oscillations that do not have a 
unique wavelength and Sk does not show a sharp cut-off. 
However, the properties are qualitatively the same as in 
the translational-invariant case. Both quantities vanish 
when Jp is reached, owing to the fact that there are no 
density fluctuations within the plateau. The character- 
istic length scales become larger as the Mott plateau is 
approached. Close to Jp we observe clear indications for 
the developing Mott plateau ^3] • The correlations of the 
density fluctuations are suppressed around the center of 
the trap leading to a local minimum of Co, r at r = 0. 
This is accompanied by a decrease of the slope of Sk ■ 

Conclusion. For the ID strongly interacting Bose gas 
in an optical lattice we have identified characteristic os- 
cillations of the density-density correlation function with 
length A. This can be used as a measure for the distance 
of the system from the MI state: the length A diverges in 
units of the lattice spacing with the density n as 1/(1 — n) 
when we approach the MI. This phenomenon is related to 
the behavior of the static structure factor Sk with charac- 
teristic wavevector k* = 2n/X. Sk is linear and vanishes 
for |fc| > k* . k* itself vanishes continuously as the MI is 
approached and Sk — in the MI phase. This behavior 
also survives qualitatively if a weak parabolic potential 
is applied to the interacting Bose gas. In particular, the 
static structure factor is strongly suppressed if a large 
fraction of the Bose gas is in the MI state. 
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